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I.  INTRODUCTION 


!  „ 

-  t 

The  shielded  suspended-substrate  line  is  a  transmission  medium  useful  for 
radar  and  microwave  circuits  in  the  Ka  band  frequency  range,  30  -  40  GHz. 
In  order  to  utilize  this  transmission  medium  in  the  construction  of  microwave 
circuits  and  filters,  it  is  necessary  to  have  valid  circuit  models  for  typical 
discontinuities  such  as  the  series  gap  in  line,  open-ended  stub,  and  a 
discontinuous  change  in  width.  These  discontinuity  characteristics  may  be 

deduced  on  the  basis  of  a  calculation  of  their  scattering  coefficients  for 

■j’hi*  {k&s  t'S 

propagating  waves  in  the  medium.  The  approach  chosen  in  -the^resgfftrwwk, 
however,  is  that  of  placing  the  selected  discontinuity  in  an  open-ended  strip 
resonator  in  the  transmission  medium.  An  appropriate  circuit  model  is 
adopted  for  the  given  discontinuity,  and  its  circuit  elements  are  then  deduced 
on  the  basis  of  the  perturbation  of  the  resonance  frequency  of  the  strip 
resonator  which  is_induced  bv  thejmodel.  ^ 

The  boundary-value  problem  associated  with  the  microstrip  resonator 
structure  has  been  approached  in  a  rigorous  manner  based  on  a  full-wave 
analysis  which  utilizes  a  process  of  solving  the  electromagnetic  (EM)  boundary 
value  problem  with  inclusion  of  all  the  field  components.  /^k\\  analysis 
presented  by  Itoh  and  Uwano  was  employed  in  the  present  work  [Refs.  1,  3]. 

The  characteristic  equation  for  resonant  frequencies  was  obtained  by  use  of 
Galerkin's  technique  applied  in  the  spectral  or  Fourier  transform  domain.  The 
set  of  algebraic  equations  among  the  transformed  quantities  thus  produced  is 
called  the  Green's  function  relations  in  the  transform  domain. 
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II.  THEORETICAL  ANALYSIS 


A.  FIELD  EQUATIONS 

The  microstrip  resonator  to  be  analyzed  is  shown  in  Fig.  1.  A 
rectangular  strip  of  width  2w  and  length  2 1  is  placed  on  the  suspended 
substrate.  The  sides  and  the  top  of  the  structure  are  surrounded  by 
conducting  shield  walls.  Thus  the  entire  structure  is  considered  to  be  a 
suspended-strip  resonator  located  in  a  partially  filled  waveguide  with  end 
walls. 


Figure  1.  End  and  top  view  of  Suspended-Substrate  Microstrip  Line 


It  is  assumed  that  the  thickness  of  the  conducting  strip  is  negligible  and 
that  all  the  media  and  conductors  are  lossless.  For  simplicity,  the  strip  is 
assumed  to  be  symmetrically  located,  although  the  present  method  of  analysis 
can  be  easily  extended  to  more  general  cases. 
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The  fields  existing  in  the  structure  shown  in  Fig.  1  are  the  superposition 
of  TE  (to  z)  and  TM  (to  z)  fields,  which  can  be  expressed  in  terms  of  two 
types  of  scalar  potentials  <J>e(x,y,z)  and  <j>h(x,y,z),  where  the  superscripts  e  and 
h  denote  electric  and  magnetic,  respectively  [Ref.  1]. 


d%e 

Ezi(x,y,z)  =  kfV  +  — - 

dz2 

(la) 

d2<j>:h 

Hzi(x,y,z)  =  k^  + 

dz2 

(lb) 

d2ty.*  4ih 

Exi(x,y,z)  =  M  ^ 

OX.OZ  dy 

(lc) 

d%h  aj>5e 

Hxi(x,y,z)  =  + 

oxaz  ay 

(Id) 

where  i  =  1,  2,  3  designates  the  substrate  or  the  air  region,  as  shown  in 
Fig.  1. 


ki2  =  k32  =  k02  =  n0e0u? 

(le) 

k22  =  w2ur  =  nTeTk02 

(If) 

el  =  e3  =  fo 

(lg) 

c2  =  6rc0 

(lh) 

Mi  =  M3  =  Mo 

(li) 

M2  =  MrMo  =  Mot  (Mr  =  i) 

(lj) 

where  w  is  the  operating  frequency  and  e0  and  /i0  are  the  free-space 
permittivity  and  permeability. 


B.  BOUNDARY  CONDITIONS 


Now  applying  boundary  conditions  at  the  bottom  in  region  1,  tangential 

electric  field  and  normal  magnetic  field  components  must  be  zero. 

d 


Ezl(x,0,z)  =  0, 

— Hzl(x,0,z)  =  0 
dy 

d 

(2a) 

Exl(x,0,z)  =  0, 

— Hxl(x,0,z)  =  0 
dy 

(2b) 

At  the  interface  between  regions  1  and  2,  the  tangential  field  components  must 
be  continuous. 

At  y  =  ht: 


Ezi(x,hi,z)  =  Ez2(x,h„z) 

(3a) 

Exi(x,hj,z)  =  Ex2(x,h1,z) 

(3b) 

Hzi(x,h!,z)  =  Hz2(x,h1,z) 

(3c) 

Hxi(x,hj,z)  =  Hx2(x,h„z) 

(3d) 

At  the  interface  between  regions  2  and  3,  tangential  electric  field  components 
must  be  continuous. 

At  y  =  h2: 

Ez2(x>h2>2)  =  EZ3(x>h2>z)  (4a) 

Ex2(x,h2,z)  =  Ex3(x,h2,z)  (4b) 

Also,  the  electric  fields  at  y  =  h2  will  exist  only  on  the  air/dielectric 
interface,  and  can  be  expressed  as, 


fEs(x)e^, 

a  >  |  x  |  >  w, 

Ez2(x,h2,z)  = 

(4c) 

[  0  , 

otherwise 

fEv(x)e12, 

a  >  |  x  |  >  w, 

II 

'eT 

CO 

X 

w 

(4d) 

l  0  , 

otherwise 

Similarly,  tangential  magnetic  fields  at  h2  must  be  discontinuous  by  the 


corresponding  surface  current  densities. 


Hz2(x,h2,z)  -  Hz3(x,h2,z)  = 

'j^xje12,  1  x |  <w, 

0,  otherwise 

(4e) 

Hx2(x,h2,z)  -  Hx3(x,h2,z)  = 

j2(x)e')z,  I  x |  <w, 

0,  otherwise 

(4f) 

At  the  top  in  region  3.  tangential  E-field  components  must  be  zero. 

At  y  =  h3: 

Ez3(X’h3>Z)  '=  °» 

0Hz3 

(x,h3,z)  =  0 
dy 

(5a) 

Ex3(x,h3,z)  =  0, 

dHxZ 

(x,h3,z)  =  0 

dy 

(5b) 

The  final  boundary  conditions  occur  at  x  =  ±  ~  where  the  tangential  E-field 

components  must  be  zero  in  all 

regions. 

At  x  =  ±-|-: 

Ezi(Hby’2)  =  0 

(6a) 

Ey  i  (±— 2~ jYiZ)  =  0 

(6b) 

C.  SPECTRAL  DOMAIN  APPROACH. 


The  boundary  value  problem  is  solved  in  the  Fourier  transform  or  spectral 
domain.  This  leads  to  the  Green's  function  equations  associated  with  the 
structure.  This  analysis  was  first  carried  out  by  Itoh  and  Mittra  and  was 
successfully  applied  to  a  number  of  problems  in  microwave  integrated  circuit 
structures  [Ref.  7]. 

1.  Potential  Equations 

The  quantities  in  Eqs.  (1)  are  transformed  into  spectral  domain  via 
the  Fourier  transform; 


00  d>  , 

$f\n,y,0)  =  f  dzf  dx$!’h(x,y,z)ejknx#z  (7a) 

“00  “a  n  =  1,  2,  .... 

where  /?  is  the  Fourier  transform  variable,  kn  is  the  discrete  transform  variable 

defined  by, 


(n  -  -~)ir 

k„  = - 1 - for  Ez  even, 


-Hz  odd(in  x)  modes. 


a 

n«7r 

kn  = - for  Ez  odd,  -Hz  even(in  x)  modes. 

a 

The  transforms  of  field  quantities  are  now, 


Ezi(n,y,/?)  =  kc| 

Hzi(n,y,/?)  =  kc| 

^  h 

Exi(n,y,/?)  =  -kn4ie  -  joyq— — 

4ie 

Hxi(n,y,/?)  =  -kn/3fih  +  jwcj— — 

Kl  =  ki2  +  r2,  (r  =  -jA  ki2  =  ^/v5) 

kcf  =  kj2  -  01  =  -  0* 

where  i  =  1,  2,  3,  for  each  of  the  three  regions  defined  in  Fig.  1. 


(7b) 

(7c) 

(8a) 

(8b) 

(8c) 

(8d) 

(Se) 

(80 
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(9a) 


The  transforms  of  scalar  potentials  satisfy, 


d2d)f’h 

— - 7.2$e,h  =  o 

dy2 


where  7j2  =  kn2  -  kc2 


7i2  =  732  =  K2  -  kcl,3  =  kn2  +  -  k02 

7 22  =  kn2  -  kc2  =  kn2  +  P  "  «  r^rk02 


The  solutions  for  these  homogeneous  differential  equations  are  well  known 
and  can  be  described  in  the  general  form  as  below, 

<j>i  =  Cjcosh^  +  C2sinh,)y.  (9b) 


Now  the  boundary  conditions  (2),  and  (5)  are  applied.  The  cosh  term 
vanishes  in  electric  field  and  the  sinh  term  vanishes  in  magnetic  field  in 
regions  1  and  3.  The  potentials  are:  (see  details  in  Appendix  A), 


Region  (1) 


<J>je  =  Asinl^y 

(10a) 

$!h  =  Bcosh7jy 

(10b) 

Region  (2) 

$2e  =  Csinh72y  +  Dcosli72y 

(10c) 

(j)2h  =  Esinh72y  +  Fcosli72y 

(lOd) 

Region  (3) 

fee  =  Gsinh73(y  -  h3) 

(lOe) 

<f3h  =  Hcosh73(y  -  h3) 

(lOf) 

where 

7i2  =  kn2  "  kci  =  kn2  +  P  ~  ki2 

(10g) 

even 

Cv 

k«  =  (n  -4-)* 

- ~  - - ,  $h  odd 

(lOh) 
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2.  Green's  Function 

In  conventional  space  domain  analysis,  the  structure  may  be  analyzed 
by  first  formulating  the  following  coupled  homogeneous  integral  equations. 
The  equations  must  then  be  solved  for  the  unknown  propagation  constant  p. 

f[  Zzz(x-x\y)Jz(x')  +  Zzx(x-x\y)Jx(x')]dx'  =  Ez(x)  (11a) 

/[Zxz(x-x\y)Jz(x')  +  Zxx(x-x',y)Jx(x')]dx'  =  Ex(x)  (lib) 

where  Ez  and  Ex  are  unknown  electric  fields  on  the  boundary  at  y  =  h2,  Jz 

and  Jx  are  current  components  on  the  strip  at  y  =  h2,  and  the  Green's 

functions  Z:  Aj  =  z,x)  are  functions  of  p. 

V 

The  integration  is  over  the  strip  where  Ez(x)  and  Ex(x)  are  zero,  as  the 

strip  is  assumed  to  be  perfectly  conducting.  In  the  spectral  domain 

formulation,  however,  the  following  algebraic  equations,  instead  of  the  coupled 
integral  equations,  are  obtained.  These  equations  are  the  Fourier  transform 
equivalent  of  the  coupled  integral  equations. 

ZzAMJA)  +  Z2x(k„,h2)3x(kn)  =  E2(kn,h2)  (12a) 

2»(k„,h2)32(k„)  +  Zxx(k„,h2)jx(k„)  =  Ex(k„,h2)  (12b) 

where  quantities  with  tildes  (~)  are  Fourier  transforms  of  corresponding 
quantities.  The  Fourier  transform  is  defined  as  in  Eq.  (7). 

The  right-hand  side  of  equations  (12)  is  not  zero  everywhere  because  the 
Fourier  transform  requires  integration  over  all  x,  not  only  over  the  strip.  The 
equations  contain  four  unknowns  Jz,  Jx,  Ez,  and  Ex,  with  unknown  P. 
However,  Ez  and  Ex  will  be  eliminated  in  the  solution  process  based  on  the 
Galerkin  procedure. 
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In  the  spectral  domain,  the  boundary  conditions  (3),  (4)  are  now  given  by 
the  following  equations, 
at  y  =  ht 


Ezi 

= 

ez2 

Exi 

= 

EX2 

Hzi 

== 

Hz2 

Hxl 

— 

Hx2 

at  y  = 

h2, 

EZ2 

= 

Ejs3 

Ex2 

= 

EX3 

hz2 

- 

hx2 

- 

hx3  = 

where  Jz  and  Jx  are  Fourier  transforms  of  unknown  current  components  Jz(x) 
and  Jx(x)  on  the  strip  at  y  =  h2. 
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When  the  fields  are  expressed  in  terms  of  the  potentials,  these  can  be  put 
into  the  form:  [Ref.  3] 


Zzz  Zzx  h 
Zxz  ZxxJ  LJx. 


where 


Zzz  =  " 


kn2+/? 


+  CZh) 


Zxz  “  Zzx  — 


K2+(P 


-  \\ 


Zxx  - 


K2+P 


-\K%  +  ^Zhl 


7y2Ct1  +  7ylCt2 


e  CtjCt2  +  Ct1Ctg-J^+  GtoCt,-^-  +  -&i- 

1  1  1  3  7y3  2  3  7y3  7y2 

-  7z2Ct2  +  7ziCtj 

"h  ""  - : - 

7zi7z2GtiCt2  +  7zi7Z3Gt!Ct3  +  7z27Z3Gt2Ct3  +  7z2 

For  7}  >  0 

Ctj  =  cotl^d,  Ct2  =  coth72t,  Ct3  =  coth73h, 

For  7j  <  0 

Ctj  =  -ja^d,  Ct2  =  -jcot72t,  Ct3  =  -jcot73h, 

7i  =  >/  kn2  +  0*  - 

7i  7i 

7yi  .  »  7zj  —  > 

J«*i  J^i 

The  quantities  Zzz,  Zzx,  Zxz  and  Zxx  are  actually  the  Fourier  transforms  of 
dyadic  Green's  function  components. 


12 


3.  Characteristic  Equations  in  Matrix  Form 

Two  equations  in  (13)  contain  four  unknowns  Ez,  Ex,  Jz  and  Jx. 
However,  the  first  two  unknowns  Ez  and  Ex  can  be  eliminated  by  applying 
Galerkin's  method  in  the  spectral  domain.  The  first  step  is  to  expand  the 
unknown  Jz  and  Jx  in  terms  of  assumed  basis  functions  Jzm  and  Jxm  with 
unknown  coefficients  cm  and  dm. 

Jz(n,/?)  =  E c  J  (14a) 

m=i 

=.|VxmW)  (Wb) 

The  basis  functions  Jzm  and  Jxm  must  be  chosen  to  be  the  Fourier  transforms 
of  space-domain  functions  Jzm(x,z)  and  Jxm(x,z)  which  are  physically  realistic, 
and  which  are  zero  except  for  the  region  |x|  <  w  and  |z|  <  L  Now, 
substituting  (14)  into  (13)  and  taking  inner  products  of  the  resulting  equations 
with  the  basis  function  Jzi  and  Jxi  for  different  values  of  i  yields  the  matrix 
equation, 


N 

E  c„J„„ 
n>l  "  “ 


+  J„;Z, 


M 

S.djvlJdj0  =  0, 


zi^zx  "i  nr  xm 
m=i 


i  =  1,  2, 


M. 


(15a) 


/PxiZxz  j^zm  +  JxiZxx  J^xJ^  =  (15b) 

i  =  1,  2, . ,  N. 


The  right-hand  sides  of  (15)  are  zero  by  virtue  of  Parseval's  theorem,  because 
the  currents  Jzi(x),  Jxi(x)  and  the  field  components  Ez(x,h2),  Ex(x,h2)  vanish 
in  complementary  regions  of  x.  For  example,  if  the  inner  product  of  Ez  on 
the  left-hand  side  of  (13)  with  JZj(kn)  is  taken,  one  obtains 

/  Jzi(kn)Ez(k„)dkn  =  2itf  JZj(x)Ezi(x)dx  =  0 

—oo  -oo 

In  the  above,  Jzi(x)  is  zero  outside  the  strip,  and  Ez(x)  is  zero  on  the  strip. 
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Therefore,  the  final  boundary  condition  is  now  used. 
Equations  (15)  will  be  expressed  in  matrix  form  as  follows, 


N  M 


EK^ic.  + 

EK<i’2'dm 

=  0, 

...,  N. 

(16a) 

1  =  1 

m=l 

i 

=  1,  2,  •• 

N 

M 

£K(2,1)c  + 

;  <  lm  m  ' 

EKii^d. 

=  o, 

M. 

(16b) 

i=l 

m=l 

i 

=  1,  2, 

where  from  the  definition  of  the  inner  products  associated  with  the  Fourier 
transform  defined  by  (7),  the  matrix  elements  are 


Ki‘>»(k0)  =  (17a) 

00  00 

K',i’2,(k0)  =  nE/o  UtMJjpJW  (17b) 

K‘2’‘>(k0)  =  ji/o\(n,«Zj2„(n,/J)d,J  (17c) 

K!S’2|«  =  Zjj&fiW  (17d) 


A  homogeneous  system  of  equations  is  thus  obtained  in  terms  of  the 
unknown  coefficients  cm  and  dm.  In  order  that  cm  and  dm  have  nontrivial 
solutions,  the  determinant  of  the  matrix  must  be  zero,  and  hence  the  resonant 
frequency  u  is  determined  for  the  resonator  and  discontinuity  assumed. 

Equations  (16)  are  now  solved  for  the  wave  number  k0  by  setting  the 
determinant  of  the  coefficient  matrix  equal  to  zero  and  by  seeking  the  root  of 
the  resulting  characteristic  equation.  The  resonance  frequency  of  the 

suspended-stripline  resonator  is  derived  from  the  obtained  value  of  k0. 

The  accuracy  of  the  solution  can  be  systematically  improved  by  increasing 
the  number  of  basis  functions  (M+N)  and  by  solving  larger  size  matrix 
equations.  However,  if  the  first  few  basis  functions  are  chosen  so  as  to 
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approximate  the  actual  unknown  current  distribution  reasonably  well,  the 
necessary  size  of  the  matrix  can  be  held  small  for  a  given  accuracy  of  the 
solution,  resulting  in  numerical  efficiency.  Hence  the  choice  of  basis  functions 
is  important  from  the  numerical  point  of  view. 


(20a) 


Jzi(z)  =  b 


-  +  “ff)’  -P2  <  2  <  -g.  A  =  P2  -  g 

l  24 sin(“7f - 5^“)>  g<z<pI?4  =  Pi-g 

With  a  gap  present  assume  Jx  =  0,  and  J2  =  ctJzl  +  CjJ^,  with  distributions 
as  shown  in  Fig.  3.  The  transformed  distributions  are: 

Jzl  =  ^Zl(^)'^zx(^n)  (21a) 

^z2  =  4  jzx(^n)  (21b) 

Jzx(kn)  is  same  as  in  Fig.  2. 

#Pi  -  e-'rfh  _  2jsin(/?g) 


izM  =  -j - — - - I  (22a) 

ftPi  +  P2  -  2g) 

'urn  = 

4  *■  7T  +  fid  J  7T  -  /?p2 

. ?-i-§rei(^  +  /))Pl-ej(^  +  /?)g . 

7T  +  Pd  j  7T  -  /?d  j 

From  Eq.  (15a)  and  with  i  =  1, 

^zi^zz^m^zm  —  ^zl(n>^)^zz^l^zl  *k  ^zl(n>/^)^zz<'2'^z2 

and,  with  i  =  2,  this  takes  the  form 

^zi^zz^m^zm  —  ^z2(n>^)^zz(h^zl  *k  ^z2(ni^)^zzc2^z2 

These  equations  can  be  represented  in  the  matrix  form: 

rc, 


j  (22b) 


(23a) 


(23b) 


Kll  K12 
K21  K22 


LC2J 


(24) 


where, 

Kll  = 

I<21  =  E/Jz2(n,/J)Zjzl(n,/W 


K12  =  s/  Jzi(n,/?)ZzzJz2(n,/?)d/? 
I<22  =  Zf3Jtt,P)iJl2(n,M3 
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B.  BOUNDARY  VALUES 

The  actual  computation  was  carried  out  using  the  shield  and  line 
dimensions  shown  in  column  i,  in  Table  1,  which  is  typical  data  from  a 
suspended  stripline  filter  problem  initiated  by  the  NOSC  (Naval  Ocean 
Systems  Center)  Microwave  Laboratory  in  San  Diego.  This  configuration  is 
used  in  the  frequency  range  of  30  -  40  GHz. 

The  data  in  columns  2  and  3  were  used  in  analysis  of  the  microstrip  line 
with  a  gap,  to  compare  with  published  data  [Ref  5]. 

Data  2  uses  frequency  range  of  85  -  102  GHz 

Data  3  uses  frequency  range  of  31  -  40  GHz 


TABLE  1.  CONFIGURATION  DATA 

(dimensions  :  mm) 


Data 

1 

2 

3 

a 

3.2 

1.27 

3.56 

d 

0.6604 

0.318 

0.763 

t 

0.254 

0.127 

0.127 

h 

0.6604 

0.19 

0.89 

w 

0.7112 

0.6 

1.5 

g 

0.26452 

0.1  -  0.6 

0.2  -  1.6 

Er 

2.2 

2.22 

2.22 

Symbols  for  Fig.  1 
.  a  =  Shield  width 
d  =  Height  of  lower  air  layer 
t  =  Height  of  dielectric  substrate  layer 
h  =  Height  of  upper  air  layer 
w  =  Microstrip  line  width 
g  =  Gap  in  microstrip  line(center) 

Er  =  dielectric  constant  of  substrate  layer. 
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CASE2 


Maximum  beta  value 


Figure  7.  Solution  frequencies  for  the  determinant 


26 


Frequency  (Hz) 


(zh)  Xouanbajj  (zh)  fouanbaij 


M 

i— < 

On 

00 

r- 

VO 

es 

rH 

C\ 

00 

t- 

VO 

w-> 

r> 

r> 

r> 

r> 

r> 

fO 

cn 

tn 

to 

r> 

r> 

r> 

r> 

ri 

ri 

ri 

ri 

ri 

ri 

ri 

ri 

to 

to 

ri 

Figure  8.  Result  of  double  summation  I 

.  i 
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Figure  9.  Resonator  and  gap  solution  with  double  summation 


28 


Maximum  beta  values 


4.  Analysis  of  Resonator  with  Gap  Present 

a.  Resonance  stub  lengths  at  given  frequency  pairs 

The  moment  method  calculation  (Appendix  D.  part  A)  was 
carried  out  with  various  assumed  stub  lengths.  The  results  are  shown  in 
Table  2  and  Table  3. 


TABLE  2.  RESONATOR  STUB  LENGTHS  WITH  DATA  2  (SEC.  1IHB) 

Units 

length  :  mm 
frequency  :  GHz 


g/2  \f 

86(85) 

94(93) 

102(101) 

0.05 

1.211647 

(1.230606) 

1.075264 

(1.090897) 

0.959969 

(0.973572) 

0.1 

1.258331 

(1.278259) 

1.121864 

(1.137848) 

1.005066 

(1.014435) 

0.15 

1.285843 

(1.304749) 

1.147675 

(1.163574) 

1.0306S6 

(1.044349) 

0.2 

1.300917 

(1.319977) 

1.163193 

(1.179146) 

1.045738 

(1.059464) 

0.25 

1.310349 

(1.329392) 

1.172264 

(1.188179) 

1.055008 

(1.068817) 

0.3 

1.316086 

(1.334999) 

1.178024 

(1.193985) 

1.060890 

(1.074570) 

c.  Effective  dielectric  constant 

The  calculation  of  the  effective  dielectric  constant  used  the  /?  values 

which  were  computed  above  (Table  4). 
u) 


d.  Characteristic  impedance  Z0  [Ref.  6] 

The  computation  of  characteristic  impedance  uses  an  empirical 
calculation  for  the  quasistatic  value,  Z0(0)  [Ref.  5]  which  has  been  confirmed 
to  be  accurate  [Ref.  8].  This  value  is  then  corrected  to  the  wanted  frequency 


f  by  use  of  Geff(f)  in  the  expression: 

z0(o) 

Z0(f)  = -  (Appendix  D.  part  B.  data  file) 

V  eeff 

(1)  for  0  <  w  <  ~~  ; 

Z,  =  60  [v  +  R- <»{<>-£-  +  y  1  +  4{-~)2J 


Zo  = 


V  =  -1.7866  -  0.2035(-]}-)  +  0.4750(-|-) 

R  =  1.0835  +  0.1007(-y-)  -  0.09457(-^-) 

(2)  when  <  w  <  a  ; 

Z1  =  120tt[V  +  R{-g-  +  1.3930  +  0.6670-Ai(-g-  +  1.444)}"1] 

z  -  Zl 

Lo  -  ~zzz zr 


V  =  -0.6301  -  0.07082(-£-)  +  0.2470(-|-) 
R  =  1.9492  +  0.1 553  (-jj-)  -  0.5123(-|-) 
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Stub  length  vs.  Gap  (Data3)  Stub  length  vs.  Freq  (Data3) 


(uitu)  ipSuaj  qnj§ 

00 


(tutu)  qiSuai  qms 


Figure  10.  Stub  lengths  related  to  gap  size  and  frequency 


33 


Gap/2  (mm)  Freq  (GHz) 


IV.  DERIVATION  OF  THE  PI-CIRCUIT 
REPRESENTATION  OF  THE  DISCONTINUITY 

A.  EQUIVALENT  NETWORK 

The  gap  discontinuity  is  represented  by  a  capacitive  Pl-network.  This 
circuit  is  terminated  on  each  side  by  an  open-ended  stub.  The  Pi-circuit 
components  are  evaluated  in  terms  of  the  admittance  matrix,  [Yy],  of  the 
two-port  discontinuity,  as  shown  below. 


Figure  11.  Pi-Network 

The  admittance  of  an  open  line  segment  is: 

Y>n  =  Ya,b  =  iYotan/?4  =  strip  line  length) 
From  the  2-port  admittance  matrix, 

V=  YUV,  +  Y12V2 
I2  =  Y21Vj  +  Y22V2 
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But  at  resonance 


V,  1  V,  1 


Yu  +  Y4  1  [y22  +  Yb1  -  YJ  =  0  (25) 

*  >  « 

This  equation  is  used  to  calculate  the  equivalent  circuit  capacitances.  The 


admittance-matrix  values  for  the  equivalent  circuit  parameters  are  as  shown  in 


B.  SUSPENDED  STRIPLINE  GAP  CAPACITANCE 

Figure  13  is  the  equivalent  network  of  Fig.  1  which  shows  the 
configuration  of  the  gap  discontinuity  in  a  suspended  stripline.  This  network 
is  enclosed  in  the  central  box  in  Fig.  14.  In  the  computation  procedure,  stub 
lengths  lx  are  assumed  and  the  corresponding  frequency  of  resonance  is  found. 
Then  a  slightly  different  length  4  is  employed  to  find  the  new  resonance 
frequency.  The  resulting  two  solution  values  provide  the  necessary  data  for 
the  calculation  of  the  discontinuity  capacitances,  Cgap  and  Cpar. 


Figure  14.  Analysis  procedure 


(Y„  +  Ya)(Yu  +  Yg)  =  Y,| 
(Y„  +  y')(Y„  +  Yg)  =  Y,2 
Y,T  +  Y„(YA  +  Yb)  +  YAY„  = 


+  Y„(y'  +  yj  +  Y  Y 


A  B 


(26a) 

(26b) 
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YA(f1)  +  YA(f2) 


(27a) 


Yu  = 


y'y'  -  YY 

A  B  A  B 

T 


Y  +  Y  -  Y  -  Y 

A  B  A  B 


Y12  =  J  (Y„  +  Y  )(YU  +  Y„)  =  *  {Yu  +  Y.(fj)} 


(27b) 


TABLE  5.  GAP  CAPACITANCES  WITH  DATA  2  (SEC.  III-B)* 


g/2  \f 

86 

94 

102 

0.05 

0.1896 

0.1893 

0.2249 

0.1 

0.1417 

0.1504 

0.1667 

0.15 

0.1209 

0.1257 

0.1346 

0.2 

0.1154 

0.1148 

0.1236 

0.25 

0.1090 

0.1073 

0.1190 

0.3 

0.1017 

0.1050 

0.1101 

TABLE  6.  PARALLEL  CAPACITANCES  WITH  DATA  2  (SEC.  III-B) 


g/2  \f 

86 

94 

102 

0.05 

20.3573 

21.5289 

22.9300 

0.1 

16.8361 

17.6024 

18.5599 

0.15 

15.0875 

15.7299 

16.4805 

0.2 

14.1369 

14.6720 

15.3471 

0.25 

13.5670 

14.0813 

14.6777 

0.3 

13.2336 

13.7098 

14.2771 

*  In  Tables  5  to  8,  frequencies  f  are  given  in  gigahertz,  and  the 
tabulated  capacitances  are  in  femtofarads. 
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TABLE  7.  GAP  CAPACITANCES  WITH  DATA  3  (SEC.  III-B)* 


S/2  \f 

32 

36 

40 

0.1 

1.4789 

1.5487 

0.9855 

0.2 

1.0326 

1.1390 

0.5812 

0.3 

0.9906 

0.4080 

0.4 

0.8200 

0.8371 

0.3403 

0.5 

0.7638 

0.2266 

0.6 

0.7307 

0.7 

0.7012 

0.1925 

0.8 

0.6594 

0.1688 

TABLE  8.  PARALLEL  CAPACITANCES  WITH  DATA  3  (SEC.  III-B)* 


g/2  \  f 

32 

36 

40 

0.1 

55.8182 

60.9670 

67.4308 

0.2 

45.8325 

48.8567 

53.0624 

0.3 

40.4891 

42.7917 

46.0020 

0.4 

37.3474 

39.3528 

41.8948 

0.5 

35.3625 

37.1171 

39.3973 

0.6 

34.0614 

35.6390 

37.7064 

0.7 

33.2005 

34.6583 

36.5680 

0.8 


32.6301 


34.0021 


35.8057 


t 


C.  RESULTS 

Tables  5  to  8  summarize  the  series  and  parallel  equivalent-circuit 
capacitances  of  the  gap  discontinuity  over  the  range  of  calculated  dimensions 
and  frequencies.  This  data  is  also  shown  in  graphical  form  in  Figs.  15  to  18. 
As  may  be  seen  in  these  Figures,  both  the  series  gap  capacitance  (Cgap)  and 
the  shunt  parasitic  capacitance  (Cpar)  were  found  to  decrese  monotonically 
with  gap  length,  at  all  frequencies.  The  plots  of  capacitance  vs.  frequency 
show,  however,  that  while  the  parallel  capacitance  rises  with  frquency,  for  all 
system  dimensions  computed,  the  series  gap  capacitance  Cgap  rises  to  a 
maximum  near  the  center  of  the  K-band  region,  but  Cgap  shows  a  trend 
toward  minimum  values  at  mid-w-band  frequencies. 
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Cgap  at  f=86GHz _  „ _  Cpar  at  f=86GHz 
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Figure  16.  Capacitances  at  W-band 
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Half  length  of  gap  (mm)  Half  length  of  gap  (mm) 
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Figure  18.  Capacitance  vs.  frequency  at  W-band 
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V.  CONCLUSIONS 


This  thesis  follows  the  methods  originated  by  Itoh  for  the  analysis  of  the 
suspended  stripline  resonator  [Refs.  1  and  3],  and  for  the  equivalent  circuit  of 
the  gap  discontinuity  [Ref.  4].  Some  changes  have  been  introduced  in  these 
approaches,  however.  A  fully-shielded  enclosure  was  employed,  allowing  the 
use  of  a  finite  Fourier  transform  in  two  coordinate  directions.  This  change 
reduced  the  computation  time  significantly,  while  reproducing  the  results  of  the 
integral  transform  along  the  line  axis. 

The  perturbed-resonator  technique  permitted  the  use  of  the  strip  current- 
density  distributions  suggested  by  Itoh  [Ref.  1].  These  are  known  to  give 
accurate  results  in  the  Galerkin  method.  A  calculation  was  made  employing 
the  resonator  configuration  used  in  [Ref.  5],  which  is  the  only  known  literature 
reference  on  the  suspended  stripline  gap. 

Table  5,  Table  6,  Fig.  16  and  Fig.  18  shows  the  resulting  parameters  at 
W-band.  Table  7,  Table  8,  Fig.  15  and  Fig.  17  show  the  parameters  at  Ka- 
band.  It  is  evident  that  all  these  capacitances  are  strongly  frequency- 
dependent.  Comparing  these  data  with  data  from  [Ref.  5],  we  see  that  the 
gap  capacitances  found  in  the  present  work  range  from  0.1  -  1.5  femtofarads. 
These  gap  capacitances  are  smaller  by  an  order  of  magnitude  than  those  found 
by  Rong  and  Li  [Ref.  5],  but  the  accuracy  of  the  latter  work  has  not  been 
confirmed.  When  the  gap  is  widened,  the  capacitance  changes  at  a  rate 
comparable  with  that  in  [Ref.  5).  The  parallel  capacitances  found  here  range 
from  13  -  70  femtofarads.  This  parallel  capacitance  shows  an  inverse 
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dependence  on  gap  length,  unlike  that  of  [Ref.  5].  The  reason  for  this 
discrepancy  is  not  known,  although  it  is  felt  that  the  capacitance  variation 
found  in  the  present  work  is  correct. 

The  two-frequency  calculation  for  each  capacitance  used  here  introduces 
some  error,  but  this  should  be  small  for  small  computed  frequency  differences. 
In  any  event,  the  circuit-equivalent  capacitances  found  in  this  way  should  be 
accurate  at  a  frequency  intermediate  between  the  two  frequencies  of  the  pair 
employed  in  each  case. 

Recommendation:  In  attempting  to  use  Itoh's  method  of  perturbed 
resonators  having  the  same  frequency,  the  determinant  of  Eqs.  (17)  did  not  go 
to  zero.  The  reason  for  this  failure  should  be  determined,  so  that  an  ideally 
accurate  computation  of  equivalent-circuit  parameters  at  a  well-defined  single 
frequency  can  be  made. 
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APPENDIX  B 


FOURIER  TRANSFORMS  OF  THE  ASSUMED  CURRENT  DISTRIBUTIONS 


The  evaluation  of  the  Fourier-transformed  forms  of  the  strip  current 
densities  is  shown  in  detail  below: 

A.  CONTINUOUS  STRIP 

The  coordinate  form  of  the  current  distributions  in  Fig.  2. 

1 


Ji(x)  = 


2w  L 


1  + 


w 


1 


J2(z)  —  — —  cos— ^  , 


J3(x)  = 

w  = 


sm- 


w 

z 


7TX 
W  ’ 


2f 


,  -W  <  X  <  w 
-t<  z  <  l 

-W  <  X  <  w 

-l  <  z  <  l 


The  Fourier  transform  is  defined  to  be: 


$(n)  =  J  (KxjeJ^dx 

"CO 


(B-la) 
(B-lb) 
(B-lc) 
(B— Id) 


(B— 2) 


The  Fourier  transforms  are  calculated  as  shown: 
.00 


3,(n)  =  /  J,(x)eik»xdx. 

— 00 

J»M)  =  /*J2(z)#zdz, 
ij(n)  =  /“j3(x)ejk«xdx, 

"CO 

UP)  =  f Jatoe^dz, 


(B-3a) 

(B-3b) 

(B-3c) 


oi 


(B-3d) 


The  Fourier  transforms  of  the  basic  current  distributions  are: 
.00 


*.(”)  -  f  * r 

-oo 


J3(n)  =  / 


-oo 

00 


-oo 


w 


- 1  +  -^jeMix, 

-W  <  X  <  w, 

(B-4a) 

cos-g-^dz, 

-e  <  z  <  e, 

(B-4b) 

.'sinJS-eiMax, 

-w  <  X  <  w, 

(B-4c) 

re^dz, 

-i  <  z  <  f, 

(B-4d) 

-oo 


The  Fourier  transforms  are  calculated  as  shown  below: 
From  Eq.  (B-4a) 


Ji(»)  =  /' -w-i1  +  HHl**  =  -sr/  [* 

-oo  L  1  J  -w 

=  -ir'f/  eik»xdx  +  —J  jx|3eiknxdxl 


(B-5a) 


The  first  term  of  Eq.  (B-5a)  is  evaluated  as  follow: 
fVu'dx  =  eiknx 

-w  jkn 


w 

x=-w 


Jkn 

1 


_(ejknw 


-K™) 


= - (cosknw  +  jsinknw  -  cosknw  +  jsinknw) 

Jkn 


=  -r- sink„w. 


The  second  term  of  Eq.  (B-5a)  is  evaluated  as  follow: 


/  |x|3ejknXdx 


(B-5b) 


(B-5c) 


-w 


Use  Judv  =  uv  -  Jvdu 
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For  positive  in  Eq.  (B-5c)  we  have: 
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For  negative  x  in  Eq.  (B-5c)  we  have: 
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Now  we  add  up  together  Eq.  (B-5b)  Eq.  (B-5d)  and  Eq.  (B-5e)  into  Eq.  (B-5a) 
1  r  2sink„w  1  ,  w3ejk“w  3w2ejk«w 


5i(n)  =  — 

2w  L  k 

6wejk«w  6e 


+ 


~{  -i - 

w3  1 


+ 


w3e-jk„W  3w2e-jkn 


w 


+  j- 


6we 


k  3 

Kn 

rjknw 


k  4 

Kn 


+ 


k  4 

Kn 


+  j" 


+ 


6e 


rjknw 


-  J- 


k  3 

Kn 

1  r  2sinknw 


k  4 

Kn 


2w 


+ 


}] 


V{ 


w3eik»w  w3e“jknw  3w2e^k»'v 
- +  j - +  - 


W 


k  2 
Kn 


3w2e_jknw  6wejk»w  6we"jk«w  6ejk»w  6e'jk«w  12 


+ 


+  J- 


l.  3 


-J- 


k  3 
Kn 


k  4 
Kn 


k  4 


■  + 


k  4 
Kn 


(B-of) 


Now  we  consider  Euler's  equation 
eik»w  -  e-jk»w  =  2jsink„w 
eik»w  +  e->k»w  =  2cosk„w 
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Finally  the  current  distributions  are  found  to  be: 
From  Eq.  (B-5f)  and  Eq.  (B-6) 
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From  Eq.  (B-7)  and  Eq.  (B-8) 

j*(M)  =  3,(n)-J4M> 

2ffsinknw  ,  r  cos  fit  sin  PC 

~  J  7T2  -  (k„w)2  l  JL  pt  (pi) 2 


7T2  -  (k„w)2  1  L  /?/  (/?. 

2fsinknw  r  cos/?/  sin/ft! 

ir  -  (knw)2  t  pt  (pi) 2 


(B— 10) 


58 


B.  RESONATOR  WITH  GAP 


Asume  Jx  =  0 

Jz  =  CiJzi  +  c2Jz2 

Jzl  =  ^zl(zWzx(x) 

*^z2  =  ^z2(z)*^zx(x) 

In  both  open-ended  line  segments  of  same  length,  current  Jzl(z)  is  chosen 
to  have  the  form  shown  in  Fig.  19. 


r - .b 

| - - - 1 

-P2  “g 

— I - 1 -  z 

g  Pi 

Figure  19.  Current  distribution  of  Jzl 

Jzi(z)  =  b>  lxl  ^  a>  lzl  ^  * 

Normalization  of  Jzl(z)  leads  to: 

(Pi  ”  g)  +  (P2  ~  g)b  =  1 
1 

b  =  - 

Pi  +  p2  -  2g 

Jzi(z)  =  b  = - —— 

Pi  +  P2  “  2§  • 

Using  Jz2(z)  from  Fig.  3,  we  have: 


-  p2  <  z  <  -  g 


g  <  z  <  Pi 


(B— 11a) 


(B-llb) 


The  Fourier  transform  from  Eq.  (B-lla)  is: 


o  o  o  o  o  o  o 


APPENDIX  C. 
PROGRAM  LISTINGS 


A.  RESONATOR  FREQUENCY 
program  SSMSL 

This  program  calculates  resonator  frequency 

integer  i,  n,  m,  compare,  count 
real*8  a,  w,  1,  d,  t,  h,  Er,  f 
resal*8  kn,  beta,  omega,  rl 
real*8  pi,  EO,  MuO,  c 
real*8  gamls,  gam2s 

real*8  Jzin,  Jzkn,  Jzbe,  Jzl,  Jxkn,  Jxbe,  Jxl 
real*8  cmp,  limit,  odlf,  newf,  limf 

logical  llog 

complex*16  gamal,  gama2,  gama3,  Ctl,  Ct2,  Ct3 
complex*16  Zen,  Zenr,  Zed,  Zer,  Zhn,  Zhd,  Zhnr,  Zhr,  Zden 
complex*16  Ze,  Zh,  Zzz,  Zzx,  Zxz,  Zxx 
complex*16  Kll,  K12,  K22,  j,  skll,  skl2,  sk22,  det 
character*8  fname 
character*  14  sname 

open(l,file='data.in') 

print*,  'Enter  output  data  file  name  fn' 
read*,  fname 
open(2,  file=fnamel 
sname  =  fname  //  '.pit' 
open(3,  file=sname) 
writef2,*]  '  file  name  =  ', fname 

write(3,*j  '  file  name  =  ', sname 

parameter 
c  =  3.e8 
j  =  cmplx(0.,l.l 
pi  =  4.*atan(l.) 

EO  =  l.e— 9/(36.*pi) 

MuO  =  4.*pi*l.e-7 
limit  =  l.e-2 
compare  =  1 
count  =  1 
llog  =  .true. 
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o  o  o  o  o  o 


call  READF(w,l,a,d,t,h,Er) 
c 

write(2,22)  w,l,a,d,t,h,Er 
22  format(t5,'Unit  :  milimiter' 

&  ,/,t5, 'Micro  strip  line  width  and  length  :  ',2(el2.6,3x) 

&  ,/,t5,'Rasonator  width  :  \el2.6 

&  ,/,t5,’heights  of  layers  1,  2  and  3  :  ',3(el2.6,3x) 

&  ,/,t5,' 'Permittivity  of  substrate  Er  :  \f7.3) 

Decide  propagation  constant  Beta 
Applying  a  root  -  seeking  process 


Calculate  the  summation 


c 


c 


c 

c 

1 


c 

11 

c 

c 


c 


c 


c 

c 

c 


print*,  'Enter  the  number  of  outer  summation  loop' 
read*,  n 
n  =  20 

write(2,*l  'The  number  of  outer  summation  loop  =  ',n 

write(3,*)  'The  number  of  outer  summation  loop  =  ',n 

print*,  'Enter  the  number  of  inner  summation  loop' 
read*,  m 

write(2,*l  'The  number  of  inner  summation  iGop  =  ',m 

write(3,*j  'The  number  of  inner  summation  loop  =  ',m 

rl  =  10.*1 

initial  value 

print*,  'Enter  initial  frequency  GHz' 
read*,  f 
f  =  f*l.e9 

print*,  'Frequency  =  ',f 
write(2,*)  'Frequency  =  ',f 

omega  =  2.*pi*f 

Kll  =  cmplx(0.,0.) 

K12  =  cmplx(0.,0,j 
K22  =  cmplx(0.,0.) 

do  20  i  =  1,  n 
kn  =  (i  -  .5)*pi/a 


Jzin  =  cos(kn*w)  -  2.*sin(kn*w)/(kn*w) 

&  +  2.*(l-cos(kn*w))/(kn*w)**2 

Jzkn  =  2.*sin(kn*w)/(kn*w)  +  3.*Jzin/(kn*w)**2 

Jxkn  =  2*pi*sin(kn*w)/(pi*¥2  -  (kn*w)**2) 

Compute  summation  loop  for  beta 
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skll  =  cmplxf0.,0.) 
skl2  =  cmplx(0.,0.j 
sk22  =  cmplx(0.,0.) 
c 

do  40  im  =  l,m 
beta  =  (real(im)  -  .5)*pi/rl 

Calculate  the  gammal  and  gamma3 

gamls  =  kn**2  4-  beta**2  -  E0*Mu0*omega**2 

call  gamct(gamal, gamls, Ctl,d) 

call  gamct(gama3, gamls, Ct3,h) 

gam2s  =  kn**2  +  beta**2  -  EO*MuO*Er*omega**2 

call  gamct(gama2,gam2s,Ct2,t) 

Calculate  Ze 

Zen  ss  gama2*Ctl/Er  +  gamal*Ct2 

Zenr  =  -j/(omega*E0) 

Zed  =  Ctl*Ct2  +  Ctl*Ct3*gama2/gama3/Er 

k  +  Ct2*Ct3  +  Er*gamal/gama2 

Ze  =  Zenr*Zen/Zed 

Calculate  Zh 

Zhn  =  gamal*Ctl  +  gama2*Ct2 

Zhd  =  gamal*Ctl*gama2*Ct2  +  gamal*Ctl*gama3*Ct3 

k  +  gama2*Ct2*gama3*Ct3  +  gama2**2 

Zhnr  =  j*omega*MuO 
Zh  =  Zhnr*Zhn/Zhd 

Calculate 

Zden  =  kn**2  +  beta**2 
Zzz  =  -(Ze*beta**2  +  Zh*kn**2)/Zden 
Zzx  =  -kn*beta*(Ze  -  Zh)/Zden 
Zxz  =  Zzx 

Zxx  =  -(Ze*kn**2  +  Zh*beta**2)/Zden 
Calculate  Jzl 

Jzbe  =  4  *pi*cos(beta*l)/(pi**2  -  (2.*beta*l)**2) 

Jzl  =  Jzkn*Jzbe 
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c  Calculate  Jxl 
c 

Jxbe  =  cos(beta*l)/(beta*l)  -  sin(beta*l)/(beta*l)**2 
Jxl  =  Jxkn*Jxbe 

Calculate  Matrix  coefficient 
c 

skll  =  skll  +  Jzl*Zzz*Jzl 
skl2  =  skl2  +  Jzl*Zzx*Jxl 
sk22  =  sk22  +  Jxl*Zxx*Jxl 
c 

40  continue 

c 

Kll  =  Kll  +  skll 
K12  =  K12  +  skl2 
K22  =  K22  +  sk22 
c 

20  continue 

c 

print*,  'check  Kll  =  ',K11 
write(2,*)  'check  Kll  =  ',K11 
print*,  'check  K12  =  ',K12 
write(2,*)  'check  K12  =  ',K12 
print*,  'check  K22  =  ',I<22 
write(2,*)  'check  K22  =  ',I<22 
c 

det  =  K11*K22  -  K12*K12 
print*,  'determinent  =  ',det 
write(2,*)  'determinent  =  ',det 
c 

call  fdsort(f, det, count) 
count  =  count  -f  1 

if(llog)  then 

c  print*,  'Make  sure  opposite  sign  between  new  and  old  one' 
c  print*,  'Not  opposite  sign  then  Enter  9  ' 

c  print*,  'Opposite  sign  then  Enter  O' 

c  read*,  ijk 

if(ijk.eq.9)  then 
oldf  =  f 
go  to  1 
endif 

llog  =  .false, 
endif 
c 

cmp  =  real(det) 
limf  =  abs(f  -  oldf) 
c 

print*,  'check  limit  freq  *,limf 
write(2,*)  '  check  limit  freq  =  \limf 
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o  o 


***************************************************** 


c 


c 


c 


c 

c 

c 

c 

c 

c 

c 

c 


c 


subroutine  slval(fun,oxv,xv,newxv,sn) 
integer  sn 

rear 8  fun,oxv,xv,newxv,posxv,negxv 

if(sn  .eq.  1)  then 
if(fun  .gt.  0)  then 
posxv  =  xv 
negxv  =  oxv 
else 

negxv  =  xv 
posxv  =  oxv 
endif 

newxv  =  (posxv  4-  negxv) /2. 
sn  =  sn  +  2 
else 

if(fun  .gt.  0)  then 
posxv  =  xv 
else 

negxv  =  xv 
endif 

newxv  =  (posxv  +  negxv)/2. 
endif 

return 

end 


subroutine  fdsort(fqy,sdet,icount) 
integer  icount,  ic 

real*8  freqy(100),fqy,rtemp,  magdet,  tfqy 

complex*16  detmt(100),sdet,cteinp,  tsdet 

tfqy  =  fqy 
tsdet  =  sdet 

if(icount  .eq.  1)  then 
freqy(l)  =  tfqy 
detmt(l)  =  tsdet 
elseif(icount  .ne.  100)  then 
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do  10  i  =  1,  i count-1 
rtemp  =  freqy(i) 
ctemp  =  detmt(i) 
c 

if(rtemp  .gt.  tfqy)  then 
freqy(i)  =  tfqy 
tfqy  =  rtemp 
detmt(i)  =  tsdet 
tsdet  =  ctemp 
endif 
c 

10  continue 
c 

freqy(icount)  =  tfqy 
detmt(icount)  =  tsdet 
ic  =  icount 
c 

elseif(icount  .eq.  100)  then 
write(3,110) 

110  format(//, t5, 'Frequency', t20, 'Mag  of  det',t35,'Determinent') 

c 

do  20  j  =  l,ic 

magdet  =  sqrt(real(detmt(j))**2  aimag(detmt(j))**2) 
write(3,210)  freqy(j), magdet, detmt(j) 

210  format(t5,fl5.2,t20, 115.3, t35,2(el5.8)) 

20  continue 

endif 
c 

return 

end 
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Calculate  Zh 

Zhn  =  gamal*Ctl  4-  gama2*Ct2 
c 

Zhd  =  gamal*Ctl*gama2*Ct2  4-  gamal*Ctl*gama3*Ct3 
&  4-  gama2*Ct2*gama3*Ct3  4-  gama2**2 

Zhnr  =  j*omega*MuO 
Zh  =  Zhnr*Zhn/Zhd 
c 

train  =  -(B(m)**2*Ze  4-  kn**2*Zh)*Jzkn**2/(B(m)**2  +  kn**2) 
sum  =  trmn  4-  sum 
c 

20  continue 

c 

mag(m)  =  sqrt(real(sum)**2  +  aimag(sum)**2) 
c 

print71,  m,B(ni),mag(m),sum 
write(2,71)  m,B(ni),mag(m),sum 
71  format(2x,i2,lx,fl2.6,lx,fl2.6,lx,2({12.6,lx)) 

10  continue 

c 

do  30  m  =  2,  (max-1) 

if(mag(m).lt.mag(m-l)  .and.  mag(m).lt.magfm+l))  then 
print*,  'Suspended  stripline  BETA  =  ',B(m) 
write(2,*) 

write(2,*j  'Suspended  stripline  BETA  =  ',B(m) 
be  =  B(m) 
endif 

30  continue 

c 

Ereff  =  (be*c/omega)**2 
print*,  'Ereff  =  ', Ereff 
write(2,*)  'Ereff  =  ', Ereff 
c 

dth  =  d  +  t  +  h 
call  impd(ZO,Ereff,a,dth,t,\v) 
c 

print*,  'Z0  =  ',Z0 
write(2,*)  'Z0  =  ',Z0 
close(l) 
close(2) 
c 

write(3,95)  f,  be,  Ereff,  Z0 
95  format(2x,el2.6,2x,f8.3,2x,f9.6,2x,f9.5) 

c 

go  to  11 
c 

999  close(3) 
c 

stop 

end 
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Input  data  file  for  frequency 

'bezdata.out' 

'f31'  31 
'f32'  32 
'f35'  35 
'f36'  36 
'f39'  39 
'f40'  40 
'f41'  41 
*f85*  85 
'f86'  86 
'f93'  93 
'f94'  94 
'flOl'  101 
'fl02'  102 


Input  data  file  for  strip  line  demensions 


Data  2. 

microstip  line  width  and  length 
Waveguide  width 
the  height  of  layers  1,  2  and  3 
the  permittivity  of  substrate 

Data  3. 

microstip  line  width  and  length 
Waveguide  width 
the  height  of  layers  1,  2  and  3 
the  permittivity  of  substrate 


:0.60000  0.00000 
:1.27000 

:0.31800  0.12700 

:002.220 


: 1.50000  0.00000 

:3.56000 

:0.76300  0.12700 

:002.220 


Output  data  file 


Frequency  Beta 
.310000E+11  683.772 

.320000E+11  705.S29 

.350000E+11  772.001 

.360000E+11  794.058 

.390000E+11  860.229 

.400000E+11  887.601 

.410000E+11  909.792 

.850000E+11  1953.919 
.860000E+11  1976.906 
.930000E+11  2137.817 
.940000E+11  2160.805 
.101000E+12  2321.716 
.102000E+12  2344.703 


Ereff  Z0 

1.109129  91.81764 

1.109129  91.81764 

1.109129  91.81764 

1.109129  91.81764 

1.109129  91.81764 

1.122532  91.26783 

1.122532  91.26783 

1.204642  82.52280 

1.204642  82.52280 

1.204642  82.52280 

1.204642  82.52280 

1.204642  82.52280 

1.204642  82.52280 


0.19000 


0.89000 
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Data  file  for  input  (frequency,  beta,  Ereff,  ZO) 
'f32'  32  31  705.829  683.772  1.109129  91.81764  8 

'f36'  36  35  794.058  772.001  1.109129  91.81764  8 

*f40*  40  39  887.601  860.229  1.109129  91.81764  8 

'f86'  86  85  1976.906  1953.919  1.204642  82.5228  6 

'f94'  94  93  2160.805  2137.817  1.204642  82.5228  6 

'f02'  102  101  2344.703  2321.716  1.204642  82.5228  6 

Data  file  for  input  (resonance  stub  length) 

0.1  'gpl'  3.280426  3.423938 
0.2  'gp2'  3.425125  3.567502 
0.3  'gp3'  3.511975  3.655768 
0.4  'gp4'  3.567063  3.711716 
0.5  'gp5'  3.604022  3.748665 
0.6  'gp6'  3.629065  3.773567 
0.7  'gp7'  3.646109  3.790344 
0.8  'gp8'  3.657755  3.801590 
0.1  'gpl'  2.789856  2.900597 
0.2  'gp2'  2.929428  3.042618 
0.3  'gp3'  3.013024  3.128288 
0.4  'gp4'  3.067654  3.181847 
0.5  'gp5'  3.104645  3.218735 
0.6  'gp6'  3.129601  3.244049 
0.7  'gp7'  3.146843  3.261261 
0.8  'gp8'  3.158684  3.272965 
0.1  'gpl'  2.398560  2.490021 
0.2  'gp2'  2.529320  2.622116 
0.3  'gp3'  2.611054  2.704331 
0.4  'gp4'  2.664520  2.758667 
0.5  'gp5'  2.701723  2.794312 
0.6  'gp6'  2.726334  2.819634 
0.7  'gp7'  2.743927  2.837242 
0.8  'gp8'  2.756226  2.849289 
0.05  'g05'  1.211647  1.230606 
0.10  'glO'  1.259331  1.278259 
0.15  'gl5'  1.285843  1.304749 
0.20  'g20'  1.300919  1.319977 
0.25  'g25'  1.310349  1.329392 
0.30  'g30'  1.316086  1.334999 
0.05  'g05'  1.075264  1.090897 
0.10  'glO'  1.121864  1.137848 
0.15  'gl5'  1.147645  1.163574 
0.20  'g20'  1.163193  1.179146 
0.25  'g25'  1.172264  1.188179 
0.30  'g30'  1.178024  1.193985 
0.05  'g05'  0.959969  0.973572 
0.10  'glO'  1.005066  1.018852 
0.15  'gl5'  1.030686  1.044350 
0.20  'g20'  1.045738  1.059464 
0.25  'g25'  1.055008  1.068817 
0.30  'g30'  1.060890  1.074570 


output  data  file  for  gap  and  parallel  capacitance 
Freq(GHz)  Gap/2  CgapffF)  Cpar(fF) 
.3200E+11  .1000E-03  1.4789  55.8182 

.3200E+11  .2000E-03  1.0326  45.8325 

.3200E+11  .3000E-03  .8946  40.4891 

.3200E+11  .4000E-03  .8200  37.3474 

.3200E+11  .5000E-03  .7594  35.3625 

.3200E+11  .6000E-03  .7176  34.0614 

.3200E+11  .7000E-03  .6859  33.2005 

.3200E+11  .8000E-03  .6594  32.6301 

.3600E+11  .1000E-03  1.5487  60.9670 

.3600E-fll  .2000E-03  1.1390  48.8567 

.3600E+11  .3000E-03  .9906  42.7917 

.3600E+11  .4000E-03  .8371  39.3528 

.3600E+11  .5000E-03  .7638  37.1171 

.3600E+11  .6000E-03  .7307  35.6390 

.3600E+11  .7000E-03  .7012  34.6583 

.3600E+11  .8000E-03  .6784  34.0021 

.4000E+11  .1000E-03  .9855  67.4308 

.4000E+11  .2000E-03  .5812  53.0624 

.4000E+11  .3000E-03  .4080  46.0020 

.4000E+11  .4000E-03  .3403  41.8948 

.4000E+11  .5000E-03  .2266  39.3973 

.4000E+11  .6000E-03  .2156  37.7064 

.4000E+11  .7000E-03  .1925  36.56S0 

.4000E+11  .8000E-03  .1688  35.8057 

.8600E+11  .5000E-04  .1896  20.3573 

.8600E+11  .1000E-03  .1417  16.8361 

.8600E+11  .1500E-03  .1209  15.0875 

.8600E+11  .2000E-03  .1154  14.1369 

.8600E+11  .2500E-03  .1090  13.5670 

.8600E+11  .3000E-03  .1017  13.2336 

.9400E+11  .5000E-04  .1893  21.5289 

.9400E+11  .1000E-03  .1504  17.6024 

.9400E+11  .1500E-03  .1257  15.7299 

.9400E+11  .2000E-03  .1148  14.6790 

.9400E+11  .2500E-03  .1073  14.08  } 

.9400E+11  .3000E-03  .1050  13:7098 

.1020E+12  .5000E-04  .2249  22.9300 

.1020E+12  .1000E-03  .1667  18.5599 

.1020E+12  .1500E-03  .1346  16.4805 

.1020E+12  .2000E-03  .1236  15.3471 

.1020E+12  .2500E-03  .1190  14.6777 

.1020E+12  .3000E-03  .1101  14.2771 
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